Proximity effect in clean strong/weak/strong superconducting tri-layers 
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Recent measurements of the Josephson critical current through LSCO/LCO/LSCO thin films 
showed an unusually large proximity effect. Using the Bogoliubov-de Gennes (BdG) equations for a 
tight binding Hamiltonian we describe the proximity effect in weak links between a superconductor 
with critical temperature T c and one with critical temperature T c ', where T c > T c \ The weak 
link (N') is therefore a superconductor above its own critical temperature and the superconducting 
regions are considered to have either s-wave or d-wave symmetry. We note that the proximity effect 
is enhanced due to the presence of superconducting correlations in the weak link. The dc Josephson 
current is calculated, and we obtain a non-zero value for temperatures greater than To for sizes of 
the weak links that can be almost an order of magnitude greater than the conventional coherence 
length. Considering pockets of superconductivity in the N' layer, we show that this can lead to an 
even larger effect on the Josephson critical current by effectively shortening the weak link. 

PACS numbers: 74.45+.c,74.50.+r 



I. INTRODUCTION 

The proximity effect between a superconductor and 
a normal metal has been thoroughly investigated using 
various techniques: Ginzburg-Landau theory 2,3 , quasi- 
classical Green function methods , Gorkov equation 
methoda 5 i 6 i 7 i 8 , and tight-binding BdG methodaW 1 . 
From an experimental point of view, one of the better 
suited experiments is the measurement of the Josephson 
critical current in weak linksi^. 

In a recent experiment 1 an unusually large proxim- 
ity effect is reported, and the authors argue that it can- 
not be explained by the conventional proximity effect. 
The system used in the experiment is a c-axis oriented 
one. The c-axis Josephson critical current is measured 
through a thin film system made of doped LCO with 
T c = 25K, sandwiched between optimally doped LSCO 
with T c = 45K. The thin film is considered to be in the 
clean limit and because of the epitaxial growth of the 
films the transmission at the interfaces is close to unity 
and interface roughness is on the order of the lattice con- 
stant. In a particular setup, the LCO thin film used had 
a thickness of 100A. Fitting the critical current around T' c 
the authors extract a coherence length in the LCO film 
which is two orders of magnitude larger than expected. 
Because of this discrepancy and the observation of non- 
zero critical current for T < 30K the authors reported 
this effect to be a "giant proximity effect" . 

Although the Josephson junction has been thoroughly 
investigated in the past for both s-wave^ and d-wave 
symmetriesi^i^iSii 7 ^*^, we feel that the calculation of 
the Cooper pair leaking distance in the case of clean limit 
and superconducting weak links needs further investiga- 
tion. We are interested to observe if the leaking distance 
will be influenced by the finite critical temperature of the 
weak superconductor. 

We propose the use of the numerical solutions of the 
BdG equations in a tight binding formulation in order to 
obtain a direct calculation of the coherence length and of 



the Josephson critical current. In the clean limit the BdG 
equations are particularly easy to solve because impurity 
averaging is not required. This method is complementary 
to the quasi-classical methods used in the dirty limit, 
namely the Usadel equations22*2I. 

For coherent transport in the c-axis direction, the 
properties of the Josephson current for d-wave super- 
conductors will be similar to the properties of the cur- 
rent for s-wave superconductors. For planar interfaces, 
with z the direction perpendicular to the interfaces, the 
d-wave order parameter will have no k z dependence, 
A(fc x ,fcj,,fc z ) ~ Ao(cos(k x a) — cos(k y a)) and therefore 
will have properties similar to a superconductor with s- 
wave symmetry. When Fourier transforming the x and y 
directions, and considering an effective ID problem in the 
z direction, the d-wave order parameter will be due to an 
effective on-site interaction within each ab-plane. We will 
calculate the Josephson current in the c-axis direction for 
a 3D d-wave superconductor. We will also show calcu- 
lations of the Josephson critical current and the Cooper 
pair leaking distance for a 2D s-wave superconductor and 
for the 100 interface of a 2D d-wave superconductor. 

The giant proximity effect is observed in underdoped 
cuprates, for temperatures T > T' c for which the mid- 
dle layer is considered to be in the pseudo-gap state. 
Previous theoretical investigations of the giant proxim- 
ity effecliSS* 2 ^ considered the N' layer to be comprised 
of pockets of superconductivity. In a recent theoreti- 
cal study«2i, interstitial oxygen dopants are considered 
to modify locally the pairing interactions. The disor- 
dered dopants are enhancing the pairing interactions, 
thus increasing the size of the local gap. This was ob- 
served in recent STM experiments* 2 - in BSCCO, which 
showed that the regions of enhanced superconductivity 
are correlated with the positions of the interstitial oxy- 
gen atoms. Because of the proximity effect, the super- 
conducting pockets will be coupled and current will flow 
through percolating paths. The presence of these pock- 
ets will effectively shorten the length of the weak link and 
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the strong external superconductors will be coupled for 
values of the effective length comparable with the leak- 
ing distance. The modification of the leaking distance 
due to the finite value of T' c will have an important influ- 
ence on the effective length. Considering equally spaced 
areas of strong superconductivity with critical tempera- 
ture T c embedded in the weak superconductor with T' c 
we will calculate the critical Josephson current and find 
its dependence on the length of the weak link and on the 
volume of the embedded superconducting pockets. If one 
considers disordered regions of strong superconductivity 
in the a — b planes of a high-T c superconductor, then 
the distance between two pockets from different Cu — O 
planes will be normally distributed. The equally spaced 
pockets scenario should be the one that gives maximal 
Josephson current and will give insight about the influ- 
ence of these pockets on the current. 

This paper is organized as follows. In the next section 
we will present our method. While the BdG procedure on 
a lattice is now well known, we nonetheless include some 
details, as some "standard" approximations are included 
for clarity. The treatment of infinite surfaces will be out- 
lined. In the third section we apply the BdG equations 
to a tri-layer system and show results of the calculation 
of the order parameter, the leaking distance and the dc 
Josephson current. Both the cases of s-wave and d-wave 
symmetries of the superconducting order parameters are 
considered. We find that the proximity effect can be con- 
siderably enhanced at temperatures close to (but above) 
the critical temperature of the weak superconductor. The 
presence of randomly distributed pockets of superconduc- 
tivity in N' enhances dramatically the Josephson critical 
current and leads to a "giant proximity effect" . 

II. METHOD 

In order to describe the superconducting state we use 
the tight binding extended Hubbard Hamiltonian: 

n = - ^4^-^X14^ 

<Lij>(j id 

+ ^Uimfnii + - ^ VijUiaiij^, (1) 

i <ij>a/3 

where Uj is the nearest neighbor hopping amplitude 
which describes the kinetic energy, \x is the chemical po- 
tential used to fix the filling of the system, Ui is the 
on-site interaction, Vij is the nearest neighbor interac- 
tion and riic, — c\ a Ci a is the density operator at site i 
corresponding to spin a. 

The properties of this Hamiltonian have been studied 
previously^; it should be viewed as an effective Hamil- 
tonian with which one can describe s-wave and d-wave 
symmetries of the superconducting order parameter. For 
an s-wave superconductor we choose an attractive on-site 
interaction Ui < and no nearest-neighbor interaction 
Vij = 0, while for a d-wave superconductor we set the 
nearest neighbor interaction to be attractive Vij < and 



the on-site interaction to vanish or be repulsive. The 
interaction parameters Ui and Vij are dependent on po- 
sition, breaking translational invariance. This will allow 
us to describe interfaces between different types of mate- 
rials. 

Using the Hartree-Fock mean field decomposition this 
Hamiltonian can be transformed into a one-particle 
mean-field Hamiltonian: 

H = ]T {-tij - S i3 n)4 a c ja + ^(Ajcfjctj + h.c.) 

<ij><7 i 

+ J2^^ c lAi +c l c h) + h - c -)- (2) 

<ij> 

For planar junctions, infinite surfaces can be consid- 
ered and therefore translational invariance in the di- 
rection parallel to the surface is recovered. By doing 
a Fourier transform of the Hamiltonian in the direc- 
tion parallel to the surface we only have to solve one- 
dimensional inhomogeneous problems. For any point in 
k-space the problem becomes a one-dimensional inhomo- 
geneous problem. In the case of a 2D superconductor 
with an infinite surface along the y direction the Hamil- 
tonian becomes: 

« = EE -(i-%)4-%[M+24^M)i4 Bt cfe 

k v <ij>cr 

+ J2J2l A * + ^cosikya^fc^f + h.c. 

fey <ij> 

where i and j are now in the direction perpendicular to 
the surface and a is the lattice constant, t and Aj~, are 
the hopping amplitude and pair potential in the direc- 
tion perpendicular to the surface and i" and A|' are the 
hopping amplitude and the pair potential in the direction 
parallel to the surface. The mean-field order parameters 
are to be calculated self-consistently: 

V h 

A i' = ^E <c Kf >«w(vo> (5) 

y k 

||E« C S C ?T> + <4 C 5»' (6) 

ky 

where Aj is the s-wave order parameter, A|' is the d-wave 
order parameter of a link in the direction parallel to the 
surface and A^ is the d-wave order parameter of a link 
in the direction parallel to the surface. 

In the 3D c-axis geometry, the surface is considered to 
be in the x — y plane. After Fourier transforming in these 
directions, the Hamiltonian becomes: 
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n = J2 E -(l-^)^-%{i" + 24 l [cos(fexa)+c OS (^a)]}ct fe, ' t c* ; f ! ' 

k x k y <ij>cr 

+ + 2A * I cos ( fc - a ) - co S {k v a)]}clf^ c k t f^ + fc.a. (7) 

fcas ky 1 



The self-consistency in the order parameters is now given 
by the following equations: 



Ui 



N * N * cX 



V Ala; ky kx ky 
< c, , C- 



>, 



(8) 



A * ~ n!n v ^ < ^ 



order parameter in the superconducting regions (we do 
not), the V = solution would need only one iteration 
to converge. 

The BdG formalism allows us to calculate the dc cur- 
rent in the absence of applied voltages. In the tight- 
binding formulation the current operator is: 



> [cos{k x a) — cos{k y a)] (9) 



where i is now taken to be the site index in the z di- 
rection. Aj is the on-site s-wave order parameter while 
A[' is the d-wave order parameter which has components 
only in the x and y directions. 

We follow the standard procedure^ of introducing a 
canonical transformation of the electron operators: 



fexfey _ ST^ ik^ky , ik^ky* \ 

T — 2-^1 " '™T ^ U n Inl 

n 

Kk y _ ST^ ik x k v _ ik^ky* i 

Hi ~ 2-^i n n '«T ■ 



(10) 

(11) 



This transformation will diagonalize the Hamiltonian and 
one obtains the BdG equations for each pair of momen- 
tum vectors fc T , fc,,: 



H 



0k x ky 
i kx ky * 



I k^g ky 



-H 



)jk x ky 

k k 



k k 

yh. X h.y 

k k 



. (12) 

These equations describe the quasi-particle states in in- 
homogeneous superconductors. The BdG equations are 
equivalent to an eigenvalue problem with parameters that 
require self-consistent calculation. We start with an ini- 
tial guess for the order parameter profile and we diago- 
nalize the resulting Hamiltonian. In our infinite-surface 
setup, we need to diagonalize a one-dimensional Hamil- 
tonian for every point in momentum space. Using the 
self-consistency equations (4)-(6) we recalculate the or- 
der parameter profile. The solution is obtained when the 
difference in the order parameters between two steps is 
smaller than a desired accuracy. 

The self-consistent calculation of the order parameter 
ensures that the order parameter in the "normal metal" 
region has knowledge about the pair potential in this 
layer. If the initial guess is a step function, i.e. the order 
parameter in the middle region is zero, after one iter- 
ation the pair amplitude will become non-zero because 
of the proximity effect. In the case U' = 0, the order 
parameter will remain zero: A ~ U- < Cj-fCjj >, while 
for the case U' < the new order parameter has a fi- 
nite value throughout the layer. If we were to fix the 



2~^i ^ ( C icr C j<? C jcr C i<r)- 



(13) 



The expectation value of the current will be non-zero only 
if the order parameters in the two superconducting layers 
have different phases. For the mean-field Hamiltonian 
with s-wave order parameters one gets for the continuity 
equation: 



drii 
< ^ > 



= < [H,th] > 

= XA'(< C L C j° > - < C ]* C i* >) 

(7 

+ < c a c lT > A*- < cf^ > A, (14) 



If the order parameter is calculated self-consistently 
Ai ^< Cjj_Cif >, then we recover the continuity equation, 
< ^ >=< Jij >. Otherwise, if the order parameters 
are not calculated self-consistently but set to a desired 
value (the case of hard boundary) then the last two terms 
in Eq. (14) can be seen as current source terms. 

In our calculation the coherence factors, Uk and Vk , are 
complex numbers and they will give the magnitude and 
the phase of the order parameters. The magnitude of the 
order parameters is calculated self-consistently and after 
each iteration the phase of the external layers is set to 
a desired value. For the calculation of the dc Josephson 
current we restrict the phase of the two superconductors 
to a desired phase difference, while for the weak link, the 
phase is calculated self-consistently. 



III. RESULTS FOR SN'S 

For the SN'S tri-layers the interactions are only on- 
site attractive interactions. The value of the parameter 
Ui will set the magnitude of the order parameter through- 
out the sample. We consider the following setup (Fig. 1), 
Ui = U for < % < A and B < i < N, while Ui = U' 
for A < i < B. In this particular case Vij is vanishing, 
because we ignore the d-wave symmetry. The value of U' 
is chosen so that \U'\ < \U\, allowing us to describe the 
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N' material with a lower critical temperature T' c . For 
temperatures greater than T' c the a — b region cannot 
sustain superconductivity by itself. The order parame- 
ter will leak from the stronger superconductors, and the 
characteristic length is called the "leaking" distance. 
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FIG. 2: (Color online) The pair amplitude at L/2 as a function 
of L for different temperatures above T c for the SN'S system 
with U' = -2t and U = -41 



FIG. 1: (Color online) The pair amplitude profile through a 
Josephson junction. The pair amplitude is shown for U = — 3t 
and U' — 0. The regions 0— A and B—N are superconducting, 
while the region A — B is a superconductor above its critical 
temperature. The interfaces are considered to have perfect 
transmission and the whole system is in the clean limit. 



Fig.l shows the order parameter profile for U = — 3t, 
V = -2t and T = 0.21*. Note that T' c = 0.205i for the 
weak superconductor, while T c — 0.459£ for the strong 
one. Similar to the T' c = case the order parameter has 
an exponential dependence on distance away from the 
interface, A ~ Aoexp(—x/£,). This is true only for tem- 
peratures much larger than T' c and for distances from 
the interface greater than the coherence length of the 
stronger superconductors. The coefficient of the expo- 
nential decay is given by the leaking distance, £. In the 
normal metal case (T' c = K) the clean limit leaking 
distance is inverse proportional to the temperature: 



(15) 



For the case of a weak superconductor, the relevant 
temperature scale is T — T' c . Plotting the order parame- 
ter versus distance from the interface on a semi-log scale 
(Fig. 2) for different temperatures, we can extract the 
leaking distance. As expected from the T' c = K case, 
the leaking distance is decreasing with increasing tem- 
perature. 

If we plot the order parameter as a function of temper- 
ature for different lengths (L) of the weak link (Fig. 3), 
we can observe two main effects. First, at T — K, the 
proximity effect will modify the order parameter at L/2. 
For lengths smaller than 10 lattice constants this effect 
is important. Because the N' layer is superconducting 
at T = OK, the main length scale in this layer is the 
superconducting coherence length £ = Th^p/A. The sec- 
ond effect is observed at temperatures close to T' c . If the 



0.2 



0.1 





L=4 

L-8 






L=16 

L 32 

L=~ 









0.2 



0.4 



T/t 



FIG. 3: (Color online) The pair amplitude at L/2 as a function 
of temperature for different L for the SN'S system with U' = 
-It and U = -41 



N' layer was not connected to the superconducting lay- 
ers, then, according to the mean-field behavior, the order 
parameter would vanish at T' c . For temperatures higher 
than T' c the N' layer cannot sustain superconductivity by 
itself. It is only in the presence of the S layers, that the 
order parameter at L/2 has non-zero values. Note that 
the length L for which we obtain non-zero values of the 
order parameter above T' c is much larger than the value of 
the length beyond which effects are unobservable atT = 
K. In Fig. 4, we compare the order parameters at L/2 
for two cases: U' = and U' = —2t. We observe that 
the order parameter for the case U' = —2t is larger and 
that close to T' c the discrepancy is enhanced. This is a 
clear indication that the Cooper pair leaking distance is 
larger in the case of a non-zero T' c . 

In order to investigate further the dependence of the 
leaking distance on the magnitude of the superconduct- 
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FIG. 4: (Color online) The pair amplitude at L/2 as a function 
of relative temperature for different L for the SN'S system for 
U = -At, U' = and U = -At, U' = -It. 




1/T (units of 1/t) 

FIG. 5: (Color online) The leaking distance as a function 
of 1/T for different interaction parameters U' for the SN'S 
system with U — —At. The vertical dashed lines represent 
the inverse of the critical temperatures for the corresponding 
U' parameters: T C (U' = -1.5t) = 0.104f, T C (U' = -2t) = 
0.205i, T C {U' = -3t) = 0.46£. 



captured in calculations that are not self-consistent. 

The proximity effect can be observed either by growing 
superconducting thin films on top of normal metals and 
measuring the critical temperature of the system, or by 
forming a Josephson junction and measuring the Joseph- 
son current through the weak link. If the two supercon- 
ducting sides are not coupled then there is no Josephson 
current. As we bring the superconducting sides closer to 
one another, the proximity effect will influence the value 
of the order parameter in the TV' layer. A non-zero value 
of the order parameter throughout the whole system will 
result in a non-zero value of the Josephson current. 

The BdG equations are well suited for calculating the 
dc Josephson current. In order to have current between 
the two superconducting sides, the order parameters in 
the two sides have to have different phases. 
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FIG. 6: (a) Phase profile and (b) dc Josephson current as a 
function of position. The phase is calculated self-consistently 
only in the middle layer and the continuity equation is satis- 
fied only in this layer. 



ing correlations in the N' layer, in Fig. 5, we summarize 
the extracted leaking distance obtained for different pa- 
rameters. The T' c = OK line(dashed) is inversely propor- 
tional to the temperature, as expected. For T' c > K the 
leaking distance is diverging at T' c \ this, of course leads 
to a giant proximity effect near these temperatures, as 
the figure visually demonstrates. Another feature of the 
calculation is that for any given temperature, a higher 
T' c will result in a larger leaking distance. For repulsive 
on-site interactions, U' = +2t, the leaking distance is 
even smaller than in the normal metal case. This is a 
clear demonstration of the fact that interactions in the 
N' layer will influence the way Cooper pairs leak from 
the superconducting side. Such "feedback" will not be 



In our calculation we fix the phases of the order param- 
eter on the S layers, and our self-consistent calculation 
will give the magnitude and the phase of the order pa- 
rameter in the N' side. The results of such a calculation 
are shown in Fig. 6a and 6b. The phase of the order pa- 
rameter in N' will vary continuously from (f>i e ft to 4> r i g ht 
and the dc Josephson current will be constant throughout 
the layer. An interesting case is the one where A(f> — n, 
for which there is a phase-slip point at L/2. Right at 
the phase-slip the order parameter vanishes. In order to 
extract only the proximity effect from the current cal- 
culation, we need to find the phase difference for which 
the current is maximal. For a point contact Josephson 
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junction the current has the following behavior—: 

./ = J m sin(A(j)), (16) 

while for a long junction it deviates from the sinusoidal 
behavio»2S. 
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FIG. 7: (Color online) The dc Josephson current in the middle 
layer as a function of temperature for different lengths L of the 
weak link for the SN'S system with U = — 3* and U' = — 1.5*. 
The arrow represents the critical temperature of the middle 
layer, T C (U' = -1.5f) = 0.104*. 

We calculate the dc Josephson current for different 
lengths of the weak link and for different temperatures 
for U' = — 1.5t and U = — 3i. The results are summa- 
rized in Fig. 7. When the two superconducting layers 
(S) are close together the proximity effect modifies the 
magnitude of the order parameter at L/2 in the N' layer. 
A large value of the order parameter will give a large 
value for the current. As L increases the order parame- 
ter decreases exponentially. This results in a decay of the 
current as a function of L. The main result is that the be- 
havior of the Josephson current as a function of L and T 
reflects the existence of a leaking distance larger than the 
one expected from a normal metal. For U' = —1.5* and 
T = 0.125*, the normal metal gives a leaking distance of 
£o ~ 2a while the self-consistently calculated one gives 
a value of £ ~ 7a. This is seen in Fig. 7, where for 
L = 16a the current is non-zero for temperatures close 
to but greater than T' c , and it has a linear dependence on 
temperature near T' c . 

As shown in previous attempts to explain the "giant 
proximity effect" , the presence of pockets of supercon- 
ductivity in the N' layer will greatly enhance the current 
through the system, even for long weak links. Coupled 
with the enhancement of the leaking distance around T' c 
the presence of the superconducting pockets will effec- 
tively decrease the length of the weak link. We consider 
equally spaced superconducting areas with on-site inter- 
actions of strength U — —4* embedded in the weak link 
with interaction strength U — —2*. In Fig. 8 we show 
the Josephson current for different lengths of the weak 



link and with superconducting pockets occupying a vol- 
ume percentage p — 0.2 of the weak link. The size of the 
considered pockets is one lattice site. The effect on the 
Josephson current is drastic — the current has non-zero 
values well above T' c . We also notice that for this volume 
of embedded superconductivity, the current has a weak 
dependence on the length of the junction. 
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FIG. 8: (Color online) The dc Josephson current in the middle 
layer as a function of temperature for different lengths L of 
the weak link for the SN'S system with U = —4* and U' = 
—2t. Equally spaced areas of superconductivity in the N' 
layer are considered. The percent volume of the pockets of 
superconductivity with U = —4* is p — 0.2. 
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FIG. 9: (Color online) The dc Josephson current in the mid- 
dle layer as a function of temperature for different percent 
volumes of embedded superconductivity in N' for the SN'S 
system with U = —4* and U = —2*. The length L of the 
weak link is L — 40a. 

The strength of the coupling between the exterior su- 
perconductors will be given by the volume of these pock- 
ets. This is seen in Fig. 9, where we plot the Josephson 
current for a weak link of length L — 40a as a func- 
tion temperature for different percent volumes of strong 
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superconducting pockets embedded in the weak super- 
conductor. As expected, for p = 0.0 the junction is too 
long to couple the strong superconductors and the cur- 
rent vanishes above but very close to T' c . Increasing p, 
the junction will effectively shorten and the two exterior 
superconductors will be coupled well above T' c . 

IV. RESULTS FOR DN'D 

The d-wave symmetry of the order parameter is 
attained if we consider nearest-neighbor interactions 
Vu+8 = —Vi and vanishing or repulsive on-site interac- 
tions. The d-wave order parameter is calculated in the 
following way: 

A d (i) = i(A a (i) + A_,(i) - A y (i) - A_ y (i)), (17) 

where A x describes superconducting correlations in the 
x direction. In a similar manner, as detailed in the SN'S 
case, we set up Vi so that Vi — V for < i < A and 
B < i < N, while V$ = V for A < i < B. This will 
allow us to describe a weak link with a non-zero critical 
temperature. 
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FIG. 10: (Color online) D-wave order parameter at L/2 as a 
function L for different temperatures - 100 d-wave case with 
V = -At and V = -2t. 



For the 100 interface (between the a — b planes of a 
high-T c superconductor), the dependence of the order pa- 
rameter is very similar to the s-wave case. Fig. 10 shows 
the semi-log plot of the order parameter as a function 
of distance from the interface for different temperatures. 
Again, we can observe the exponential decay and define 
the leaking distance £. The dependence of the order pa- 
rameter on temperature for different lengths of the weak 
link is shown in Fig. 11 and the two manifestations of the 
proximity effect are seen. First at T = OK the order pa- 
rameter is modified if L/2 is of the order of the supercon- 
ducting coherence length in the N' layer. Secondly, above 
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FIG. 11: (Color online) D-wave order parameter at L/2 as 
a function of T for different lengths - 100 d-wave case with 
V = -At and V = -2t. 
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FIG. 12: (Color online) Leaking distance as a function of 
inverse temperature for different interaction strengths - 100 d- 
wave case with V — —At. The vertical dashed lines represent 
the inverse of the critical temperatures for the corresponding 
V' parameters: T C (V' = -2i) = 0.4i, T C (V' = -3t) = 0.67t. 



T' c the order parameter decays with increasing tempera- 
ture but has non-zero value even if L/2 is greater than 
the conventional leaking distance defined by the T' r = OK 
case. The self-consistently calculated leaking distance is 
shown in Fig. 12. Similar to the s-wave case it diverges 
at T' c and, for the same temperature, larger interactions 
in the weak superconductor will increase the leaking dis- 
tance. 

The coherent transport in the c-axis direction will be 
described by the hopping amplitude in the z direction, 
t 1 - = 0.5i". Fig. 13 shows the Josephson critical current 
as a function of temperature for different lengths of the 
weak superconducting layer. The general behavior is sim- 
ilar to the one of the s-wave junction. The d-wave order 
parameter has no k z dependence, and thus in the z direc- 
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FIG. 13: (Color online) The c-axis dc Josephson current in 
the middle layer as a function of temperature for different 
lengths L of the weak link for which V = —At and V = — 2t. 
The c-axis hopping amplitude is t 1 " = 0.5i". 



tion it is an effective on-site order parameter. For short 
weak links the current does not vanish abruptly above 
T' c but rather has a smooth dependence on temperature. 
This dependence on temperature above T' c shows that the 
length of the weak link is comparable to the leaking dis- 
tance in this layer. The increase of the leaking distance 
due to the finite T' c cannot explain by itself the observed 
"giant proximity effect" . It is only the conjunction with 
the presence of disordered pockets of superconductivity 
in the weak link that makes this effect possible. The cal- 
culation of the Josephson current in the presence of the 
disordered pockets from the previous section stands also 
for the c-axis geometry. The extra dimension will only 
affect the necessary volume of superconductivity needed 
to observe a "giant proximity effect" . 



the proximity effect induced by the exterior supercon- 
ducting layers in the "normal metal" interior layer. We 
observed that, in agreement with previous calculations, 
the order parameter has an exponential decay behav- 
ior, the characteristic decay length being the leaking dis- 
tance. In both s-wave and d-wave cases the leaking dis- 
tance is only dependent on the properties of the N' layer. 
For T' c = OK (normal metal) for both s-wave and d-wave 
symmetries the leaking distance is inversely proportional 
to the temperature. If T' c > OK, the leaking distance 
diverges at T' c and at the same temperature larger at- 
tractive interactions in the middle layer will increase the 
leaking distance. Essentially, the BdG formalism pro- 
vides a means for the normal layer to feel pairing fluctu- 
ations above its critical temperature, T' c . These are not 
spontaneous, in that they arise from an 'applied' pairing 
field produced by the outer layers. This accounts for the 
much higher leaking distance for a weak superconductor. 

We also calculated the dc Josephson current, and ex- 
tracted the maximum value. We observed that the cur- 
rent has a non-zero value for lengths of the weak link 
much larger than the "conventional" leaking distance, 
and for temperatures well above T' c . Although the di- 
vergence of the leaking distance at the critical tempera- 
ture of the N' layer enhances the Josephson current for 
temperatures above T' c , it is not enough to explain the ex- 
perimental measurement of the "giant proximity effect"—. 
As prompted by previous attempts to explain the "giant 
proximity effect"^2*2£, we considered areas of supercon- 
ductivity with critical temperature T c > T' c , which are 
embedded in the N' layer. Further enhancement of the 
Josephson current is observed. Depending on the vol- 
ume of the superconducting pockets, non-zero values of 
the Josephson current are obtained even for temperatures 
T > 2T' C . These results form the basis for a qualitative 
understanding of the giant proximity effect observed by 
Bozovic et alii 



V. SUMMARY 

In summary, using a tight binding formulation of the 
extended Hubbard Hamiltonian, we solve the BdG equa- 
tions for a system composed of three layers: two super- 
conducting layers (either s-wave or d-wave), and a weaker 
superconductor sandwiched in between. We examined 
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